Polaron tunneling dynamics in the DNA molecule 
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The formation of polaron and its migration in a DNA chain are studied within a semiclassical 
Peyrard-Bishop-Holstein polaron model. Comparing the energetics of the polaron system found 
from the quantum chemical and semiclassical calculations, we extract the charge-phonon coupling 
constant for poly DNA sequences. The coupling constant is found to be larger for the G-C than for 
the A-T pairs. With this coupling constant we study tunneling in the DNA molecule. The rates 
and the nature of tunneling have strong dependence on the DNA sequence. By changing the trap 
positions in the molecular bridge the tunneling rate can be varied up to seven orders of magnitude. 



The discovery of conductance in DNA has attracted 
many researchers to investigate the transport properties 
0,12 Si, IE IS 0,110 of DNA. For the (G-C)(A-T)jv(G- 
0)3 DNA sequences the mechanism of charge transfer is 
more or less clear 0, S an d is described by the com- 
petiton between tunneling and hopping transfer. But 
for the poly and mixed DNA sequences the experimental 
data observed by different groups are often contradictory. 
In some experiments a high conductivity was obtained 
, while in others the conductivity was rather low 
In the works devoted to simulation of charge trans- 
fer in the DNA molecule within the tight-binding Hamil- 
tonian [IS [ll[ or the system of kinetic equations H , no 
explanation of this phenomena was found. Within these 
models the charge transfer integral between the nearest 
base pairs and the energy gap between the states were 
the key parameters. At the same time the models did not 
take into account consistently the effect of geometry fluc- 
tuations or phonons on the charge transfer processes in 
the DNA molecule. This is despite the fact that already 
in 1956 Marcus pointed out [12J that geometry fluctua- 
tions can activate and strongly affect an electronic tran- 
sition between two states. Therefore, the polaron model 
[S IS US 13 , which takes into account all these param- 
eters should be invoked for an adequate investigation of 
the charge transfer in DNA. 

In this paper, we demonstrate the possibility to de- 
sign an artificial DNA molecule with semiconductor or 
insulating behaviors simply by placing a trap at the cor- 
rect points. The study is based on the analysis of charge 
transfer from the donor to the acceptor through a molecu- 
lar bridge composed of the potential barriers (A-T pairs) 
and the charge traps (G-C pairs). We show that in the 
mixed DNA molecules the transfer rate of the charge 
strongly depends on the sequences, i.e., on the positions 
of the traps between the donor and the acceptor. This po- 
sition determines the relation between the rate of charge 
trapping and the rate of charge escape from the trap. 
Depending on this ratio the tunneling between the donor 
and the acceptor can be described either as a sequential 
tunneling or coherent tunneling through a trap. This 
difference strongly affects the final rate of charge trans- 
fer. We show that the time of polaron tunneling can be 
changed by 10 7 times by changing the position of the 



trap. The slower transfer occurs when the polaron trap- 
ping rate is much slower than the escape rate. On the 
other hand, the fastest polaron transfer occurs when the 
trapping and the escape rates are equal. In this case the 
tunneling has the coherent nature and the polaron only 
partially occupies the trap. 

We focus here on the DNA structures where the gua- 
nine and the adenine are stacked in one DNA strand. 
In this case, there is only one degree of freedom for the 
charge transfer - the longitudinal one-dimensional charge 
tunneling through a single strand. The polaron tunnel- 
ing is described by the system of equations within the 
Peyrard-Bishop-Holstein (PBH) model [H, G3], where 
the charge motion is treated quantum-mechanically while 
the polaron tunneling classically. The Schrodinger equa- 
tion describing the dynamics of the charge within the 
DNA chain with n sites is determined by the Hamilto- 
nian 7i, which includes the tight-binding and the charge- 
lattice interaction terms. The Schrodinger equation has 
the form 



1—— 

dt 



= -Vi. 



Vi, i+ i^ i+1 +Xiyi^i-^i, (i) 



where \&j is the probability amplitude for the charge to 
be on the i-th base pair, Vi—i ^(Vt^+i) is the transfer 
integral between the nearest base pairs, Xi is the charge- 
vibrational coupling constant for the i-th site, e, is the 
on-site energy, yi determines the stretching at the z-th 
site, i.e. the displacement of the atomic structure. The 
motion of the stretching displacement yi is described by 
the Newton equation as [14| 



d 2 yi 
1 dt 2 



= -VUVi) - W'favi-x) - W'(y i+uyi ) 



1 it, 1 2 d v* 



(2) 



where m is the base pair mass, 7 is the friction param- 
eter (7 = 1 ps _1 HI])) Viaiyi) is the vibronic potential 
of the deformation energy of the hydrogen bonds in the 
base pair taking into acount the repulsive interactions of 
the phosphates and W(yi,yt-i) is the nearest- neighbor 
potential of interactions of the stacked base pairs |T3 |. 
The parameters for the deformation potential and the 
interaction potential are taken from Ref. [3| • 
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TABLE I: The values of |Ajv and Xi for the (G-C)jv and the 
(A-T)jv complexes (V<_i,i = Vi, i+ i = V<=0.1 eV). 



\Xn (eV) 



(eV/A) 



(G-C)iv = i 
(G-C)jv=2 

(G-C)at=3 

(G-C)jv = 4 
(A-T) N =i 
(A-T) w =a 
(A-TV=3 
(A-T) N =4 



0.360 
0.310 
0.265 
0.225 
0.185 
0.165 
0.140 
0.125 



1.60 
1.15 
0.90 
0.60 
1.05 
0.53 
0.40 
0.30 



"Simulation results are from Ref. [IE 



The probability amplitude and the position of 
the stretching displacement in time are evaluated 
from the self-consistent solution of the time-dependent 
Schrodinger and Newton equations. Here the 
Schrodinger equation (1) is presented as 



1+ 2H HAt 



*(t + At) 



2E HAt 



*(t) (3) 



where the Hamiltonian Tt is a function of the stretching 
displacements yi{t). The three-point difference scheme 
for the first order time derivation of the stretching pa- 
rameter has been used in Eq. ^j. For the stationary 
solution of Eqs. fl][2]), the initial occupation probability 
is assumed to be close to unity on the donor site. This 
solution is taken as the initial state (t=0) for studying 
the polaron dynamics. 

At first we analyze the equilibrium stationary pola- 
ronic states within a finite region of the DNA chain. The 
charge- vibration coupling constant Xi is the- main param- 
eter regulating the stretching of the polaron to the near- 
est sites [3] and the magnitude of y^. The value of Xi 
depends on the geometries of the DNA sites participating 
in the formation of the polaron. The shift of the state 
energy due to the polaron occupation XiUi m the absence 
of the DNA-solvent interactions can be described by the 
inner-sphere reorganization energy (« 0.5Ai in Ref. [HI) 



1 N 
-Ajv « 22xiVi, 
1 <=i 



(4) 



where N is number of sites occupied by the polaron. Re- 
cently, the exponential decrease of the inner-sphere reor- 
ganization energy with the elongation of the DNA chain 
was found within the quantum chemical calculations for 
the (G-C)jv and the (A-T)jv chains The geometry 
relaxation was found to have a maximum at the center 
of the polaron, which agrees with the results obtained 
within the PBH model [l4 |. 

The reorganization energy Atv found in Ref. [l5| are 
shown in Table I for the (G-C)jv and the (A-T)jv chains. 
From these data we can estimate the coupling constants 
Xi for different systems. The corresponding results for Xi 




FIG. 1: The wave function ^ and the lattice displacement 
yt for the polarons formed in GGGGGGGGG (solid line) and 
AAAGGGAAA (dashed line) chains. 



are shown in Table I. We have found that the coupling 
constant is smaller for the A-T base pair than for the G-C 
pair. We also determine the tendency in the dependence 
of the coupling constant on the size of the complexes; 
the coupling constant decreses with increasing sizes of 
the (G-C) at and (A-T)jy chains. 

With the values of coupling constants derived for dif- 
ferent base pairs we now study the properties of the po- 
laronic state in the poly and mixed DNA chains. In 
the poly(dG)-poly(dC) and the poly(dA)-poly(dT) DNA 
molecules, according to our calculations the polaron oc- 
cupies mostly 7-9 sites. In the mixed DNA chain, 
the polaron stretching is limited by the difference be- 
tween the coupling constants Xii and on-site energies ej, 
for A-T and G-C pairs (see Tabic [J). The results for 
GGGGGGGGG and AAAGGGAAA chains are shown 
in Fig. 1. The polaron in the AAAGGGAAA structure 
is mostly localized within the GGG due to high potential 
barriers between guanines and adenines (1.08 eV) fill ]. 

The value of the coupling constant also determines the 
polaron stretching, but its effect is strong only in the 
structure with low potential barriers. An example of 
such structures is the AAAAGAAAA chain where the 
energy gap between the A-T and the G-C is only 0.4 eV. 
The results for the AAAAGAAAA structure are shown 
in Fig. [3 (a) for two cases: (i) when the value of coupling 
constant is the same over the whole chain and is equal to 
Xi=0.6 eV/A [lj|, and (ii) when the coupling constant 
is different for G-C and A-T base pairs. Clearly, the in- 
troduction of different coupling constants Xi f° r A-T and 
G-C pairs provides stronger localization of the polaron 
within the G-C trap. 

The effect of the coupling constant on the polaron lo- 
calization in the GGGGAGGGG chain is illustrated in 
Fig. [2] (b). Actually, for this chain the polaron vibration 
mode is outside of the lattice band of the A-T site [Til ]. 
When the coupling constant Xi=0.6 eV/A is the same 
over the whole DNA chain, the vibration mode is only 
marginally delocalized. The energy of this state is —0.81 
eV, while the potential barrier between (G-C) n and A- 
T site is —0.87 eV. The polaron in this case is almost 
localized at the A-T site. When we introduce the depen- 
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FIG. 2: The wave function and the lattice displacement yi 
for polarons formed on (a) G-C for x»=0.6 eV/A (solid line) 
and xg-c = 1-6 eV/A, Xa-t = 0.4 eV/A (dashed line) and 
on (b) A-T Xi=0.6 eV/A (solid line) and xg-c=0.9 eV/A, 
Xa-t=1.05 eV/A (dashed line). 



dence of the coupling constant on the site type [Fig. [5] (b) 
(dashed line)], the energy of the state becomes —0.56 eV 
and the polaron becomes delocalized over three nearest 
G-C sites from each side of the A-T pair. 

To study polaron tunneling between the DNA traps, 
we first compare the energies of the polaronic states in 
different types of traps. Localization of the polaron in 
the (G-C)jv traps shifts the on-site energy e, to a lower 
value. This is the energy of the polaron which is the 
eigenvalue of the Hamiltonian corresponding to Eq. (|TJ) - 
This energy can also be estimated from the site energy 
Cj and the electronic energy XiVi [3 



N 



N N 

Y j e l /N=-\ N + Y.^/ N 



(5) 



For the energy difference between the polaronic states in 
different traps we have found the values —0.20 eV for G- 
C and (G-C) 2 traps and -0.43 eV for G-C and (G-C) 3 
traps, fnclusion of inner-sphere reorganization energy 
into the charge transfer model has brought down these 
values from Ae = —0.47 eV and Ae = —0.68 eV, respec- 
tively [161 ]. A direct comparison with the experimental 
results in the solvent [l7| would require evalution of the 
solvent reorganization energy [l8| . which is beyond the 
scope of this work. However, for the results that follow, 
in particular for the polaron migration dynamics, the sol- 
vent contribution perhaps is not the dominant one. 

The low energy gap between the states of the (G-C) n 
traps results in the competition between two processes 
in the mixed DNA [l9(: (i) the trapping of the polaron 
within the trap and (ii) the tunneling of the polaron be- 
tween the (G-C) m traps. To study the problem of po- 
laron tunneling between the DNA traps we have per- 
formed numerical simulation of the polaron dynamics in 
a mixed DNA chain. Here the first G-C trap is a donor 
with localized charge on it in the initial state and the 
(G-C)3 trap is an acceptor. Since the system is initially 
in the nonequlibrium state the polaron will tunnel from 
the donor to the acceptor. For the system without any 
additional traps, i.e. a DNA chain with only the donor 
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FIG. 3: The dynamics of propagation of the wave func- 
tion and the lattice displacement yi in the chains: (a) 
G(A)iGG(A) 3 GGG and (b) G(A) 2 GG(A) 2 GGG. 



and the acceptor [(G)(A)„(G)3], we have found an ex- 
ponential dependence of the tunneling rate on the tun- 
neling distance. Let t n be the tunneling time for the 
structure (G)(A) n (G)3. We then have for the normal- 
ized rate, tt/t„ = 0.6 x 10 _1 for n = 2, 1.4 x 10~ 3 for 
n = 3, and 0.6 x 10~ 5 for n = 4 (t n = 0.5 ns). These 
data are in a good agreement with the tunneling rates in 
the experimental results (20I ]. 

The new features in the polaron tunneling process is 
observed for the mixed DNA structure with additional 
trap between the donor and the acceptor. Here we study 
the system with the (G-C)2 trap in the (A-T)g molec- 
ular bridge. The dispersion of the site energies of the 
G-C base pairs within the (G-C)3 where the guanine at 
the end site has higher site energy than guanines located 
close to the sequence center [ll| , has been taken into ac- 
count. In the case of the (G-C)2 trap located close to the 
donor site (G-C) the polaron is stretched over the donor 
and the trap. As a result the polaron quickly tunnels to 
the (G-C)2 trap [Fig.|5{a)]. The polaron occupation pro- 
cess takes some time and finally the polaron tunnels to 
the acceptor. In this case the tunneling from the donor 
to the acceptor states has the sequential nature and the 
tunneling processes from the donor to the trap and from 
the trap to the acceptor are uncorrelated. 

When the (G-C)2 trap is placed exactly in the middle 
of the (A-T)6 bridge, a significant change in the polaron 
tunneling dynamics is observed [Fig.[3](b)]. In this case, 
the rate of charge tunneling from the donor to the trap 
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FIG. 4: The dependence of the occupation probability |^i| 2 
on time for the chains: (a) G(A)iGG(A) 3 GGG and (b) 
G(A) 2 GG(A) 2 GGG. 

is almost equal to the rate of tunneling from the trap 
to the acceptor. Therefore, the polaron is only partially 
localized on the trap and the final polaron tunneling from 
the donor to the acceptor is a coherent process. The 
curve for t = 10" 6 sec in Fig. 0(b) shows the occurence of 
the resonance effect due to the coincidence of the trap site 
energy with the site energy of the last guanine within the 
(G-C)3 acceptor. The rate of charge tunneling in F ig, p 
is in good agreement with the experimental results [2l| . 
where the transfer from a donor to an acceptor in similar 
systems was estimated to be 10~ 8 — 10 -6 sec. 

In Fig.[4]the dependence of the occupation probability 
of different traps within the DNA chain is shown for (a) 
G(A)iGG(A) 3 GGG and (b) G(A) 2 GG(A) 2 GGG struc- 
tures. We again see a completely different nature of tun- 
neling for different positions of the trap. When the trap 
is close to the donor, the charge transfer process is the 
sequential incoherent tunneling, i.e., when the polaron 
spends a long time within the trap. But if the trap is 
moved closer to the center of the tunneling bridge then 
the tunneling becomes coherent. However, when the po- 
sition of the (G-C)2 trap is closer to the acceptor as in 
the G(A) 3 GG(A)iGGG structure, the polaron is not lo- 
calized on the trap but tunnels directly from the donor to 
the acceptor. Since the trap does not participate in the 



charge transfer, the width of the potential barrier for the 
polaron covers the whole molecular bridge (A)3GG(A)i 
and coherent tunneling from the donor to the acceptor 
occurs in the range of t — 200 sec. This is 10 7 times 
slower than the time for coherent tunneling through the 
trap [see in Fig. 3(b) and Fig. 4(b)]. Therefore, the trans- 
fer mechanism for the G(A)3GG(A)iGGG sequence is 
similar to that for the G(A)gGGG structure. 

In conclusion, from the results of the ab initio quan- 
tum mechanical calculations, we obtained the charge- 
vibration coupling constants in the Pcyrard-Bishop- 
Holstein model for polarons formed in the (G-C)jv and 
the (A-T)^r DNA molecules. We have found that the cou- 
pling constants are larger for the (G-C)jv complex than 
for the (A-T)jv. From the calculated values of the cou- 
pling constants we have studied the energetics and the 
structure of the polaron in different DNA sequences. In 
the poly-DNA molecule, the polaron occupies nine DNA 
base pairs, while in the mixed DNA the size of the po- 
laron is strongly affected by the potential gap between 
the A-T and the G-C sites. In addition to the properties 
of the stationary polaronic state, we have also studied 
the dynamics of the polaron tunneling from the donor to 
the acceptor. We have found a very strong dependence 
of the tunneling rates on the structure of the tunnel- 
ing bridge. The position of additional traps within the 
bridge strongly affects the nature of the tunneling process 
and the rates. By changing the position we can change 
the tunneling rate upto seven orders of magnitude. To 
have the fastest tunneling rate we need to have coherent 
tunneling, i.e. tunneling to each of the traps should be 
almost equal to the escape rate from the trap. Our cal- 
culations are restricted to low temperatures. Therefore, 
the dynamics of the polaron in our study is completely 
due to quantum mechanical processes. 
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